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Abstract. 

The ASTRA Cassegrain Spectrophotometer and its automated 0.5-m tele- 
scope at Fairborn Observatory in Arizona will produ ce a large quantity of high- 
precision stellar flux distributions. A separate paper ijAdelman et al.ll2007() pre- 
sented a review of the design criteria for the system and an overview of its 
operation. This paper discusses the techniques used in the data reduction to 
final flux calibrations. 

Extraction of 1-d spectra fro m the 2-d images wil l be performed by a highly 
automated version of CCDSPEC (G ulliver fc 11 111 2002'). The characteristics of 
the CCD are automatically applied to the images, including the location of dead 
rows and hot pixels. In order to achieve the goal of better than 1% precision, 
large numbers of bias and flat field frames will be used in the reduction process. 
There will be a continual programme to monitor the image quality. Finally, 
optimally extracted spectra will be obtained, including the removal of scattered 
light and cosmic rays. 

The Earth's atmosphere has a considerable effect on the stellar flux as mea- 
sured from the surface. The principal sources of extinction, Rayleigh and aerosol 
scattering, ozone and telluric line absorption, are discussed, along with methods 
used to determine their effects on the observed spectra. Correction for telluric 
lines is the most problematical, due to their non-linear variation with airmass. 
By using a large network of constant stars to monitor atmospheric extinction it 
is possible to determine the extinction coefficients to generally better than 1% 
and to assess their temporal variability. 

The spectrophotometric observations are placed on an absolute flux scale 
by reference to stars with known values of true flux at top of Earth's atmosphere. 
These standard stars have been calibrated against terrestrial sources of known 
properties. Unfortunately, very few stars have been calibrated directly. The 
ASTRA fluxes will be calibrated against the best available Vega flux distribution. 
The constant stars used in the extinction determinations will provide the internal 
calibration network of secondary flux standards. 

The available absolute calibrations are accurate to typically 1-2%. Ulti- 
mately this uncertainty will limit the accuracy of the final fluxes of other stars. 
However, the internal precision will be significantly higher, and should more- 
accurate absolute calibrations become available the fluxes can be re-calibrated 
to higher accuracy. 
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1. Introduction 



The ASTRA Cassegrain Spectrophotometer and its automated 0.5-m telescope 
at Fairborn Observatory in Arizona will prod uce a large quantity of high-precision 
stellar flux distributions. A separate paper, ()Adelman et al.ll2007t ^. presented a 
review the design criteria for the system and an overview of its operation. In 
this paper we discuss the techniques used in the data reduction to final fiux 
calibrations. 

ASTRA will produce a vast quantity of high quality data. The reduction 
procedures will be automated as far as possible, using software sentinels to watch 
for unusual events. 



2. Reduction of 2-d Images 



The spectrophotometric observations taken by ASTRA comprise CCD frames 
containing two spectral orders. The orders are sufficiently separated that there is 
no order overlap and each order contains both the target spectrum as well as sky 
background. Spectral extraction is performed usin g a highly autom ated version 
of the CCDSPEC reduction package developed bv iGulliver fc Hilll ((2003). 

Spectral Extraction follows the procedures common to normal spectra, 
starting with correction for CCD properties, including dead rows, hot pixels. 
Next, bias corrections and flat flelding is performed using master calibration 
frames obtained from the addition of many bias and flat fleld frames. Order 
tracing is performed and PSF fitting using optimal extraction, followed by scat- 
tered light reniqval a nd the removal of cosmic rays. Details can be found in 
IGulliver fc ffilil (|2003). 

As part of the commissioning process we will characterise the extent of 
CCD fringing, determine shutter timing corrections, and assess the nature of 
any flexure of spectrograph. 

Unlike normal spectroscopy, spectrophotometry can be regarded as precise 
relative photometry at many wavelengths. Hence, the spectral extraction pro- 
cedure must recover exactly the same relative proportion of observed flux at 
every wavelength from all observed spectra. Since the spectral width depends 
on atmospheric seeing and the amount of telescope trailing, the extraction aper- 
ture must be varied to ensure the same fractional spectral energy is covered. 
The iHornd (jToSfi ) optimal extraction algorithm maintains spectrophotometric 
calibration validity. 

It is also important to preserve the observed count level and not re-bin 
onto a linear wavelength scale, since this will affect the atmospheric extinction 
determination. The output from this stage of the reduction procedure is instru- 
mental wavelengths and counts/s per wavelength bin, plus an estimate of the 
uncertainty in the counts. 

Figure ^ shows an example of a simulated extracted spectra for a 10th- 
magnitude AOV star. It shows that the 1-hour integration has produced a signal- 
to-noise (S/N) in excess of 200:1 in most of the range from AA3700-9300. 
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Figure 1. Simulated Extracted Spectrum. This is given in estimated S/N 
for a 3600-second exposure of a 10th magnitude AOV Star. 

3. Earth's Atmosphere 

The Earth's atmosphere has a considerable effect on the stellar flux as measured 
from the surface. The principal sources of extinction, Rayleigh and aerosol scat- 
tering, ozone and telluric line absorption, will be discussed, along with methods 
used to determine their effects on the observed spectra. 



Rayleigh Scattering Rayleigh scattering due to m olecules in the atmosphere is 
propo rtion al to A~^, with d etailed formulae given bv lAllenI (|l973) , iHaves fc LathamI 
(1197,^ and lFr5hhch fc Shaw (198nh . Vertical extinction is proportional to local 
atmospheric pressure. 

Aerosol Scattering Aerosol scattering is due to dust particles, salt particles, wa- 
ter dr oplets, man- made pollutants in the earth's atmosphere ()Haves Sz LathamI 
Il97,ll : iBurki et al. I ll99,l iForbes et al.lll99^ . We can model aerosol absorption 
using Angstrom's simple empirical approximation formula AqX~°'. Aerosol ex- 
tinction is quite variable, with showing diurnal and seasonal variations to a 
factor of two or more, and 0.5 ^ a ^ 1.5. 

Ozone bands The Huggins band provides the ultraviolet cut-off and the Chap- 
puis band is present in the green part of the optical spectrum, with two strong 
diffuse peaks around 5730A and 6020A. High- resolution absorption coefficients 
are available ( Burrows et al.l 119991 : IVoig3l200ll l. Ozone is concentra ted at alti- 



tudes between 10 and 35 km, and can vary significantly in a few hours ( Haves Sz Lathainl 
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1197.^ . Typically, Chappuis can contribute around 0.01 mag. of absorption per 
airmass. 

Telluric lines Several bands of discrete absorption lines, especially in the red 
and near-IR: oxygen (O2) around 7590A and 6870A, and water vapour (H2O) 
bands near 7100A, 8090A, 8920A, and 9277A. Water vapour is highly variable 
with irregular night-to-night and seasonal variations. The H2O bands do not 
follow the normal exponential absorption (Bouguer) law, but show a curve-of- 
growth effect iHaves. (,197^) . 
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Figure 2. Simulation of the typical extinction by Earth's Atmosphere, show- 
ing the relative relative contributions of Rayleigh and aerosol scattering, ozone 
and telluric line absorption. 

3.1. Extinction Law 

Extinction by the Earth's atmosphere is given by Bouguer's Law, 

mx = mxo + kxX, 

where mx is observed stellar magnitude at wavelength (A); mxo is stellar magni- 
tude above Earth's atmosphere; kx is extinction coefficient (mag. /airmass); and 
X is the airmass relative to that at the zenith. 

The total extinction coefficient at each wavelength is the sum of the various 
individual contributions: Rayleigh, Aerosol, Ozone, Oxygen and Water (see 
Fig. El). 



ASTRA Spectrophotometer: Reduction and Flux Calibrations 



5 



3.2. Air mass 



In the classical text, iHardid (|l962r i gives a polynomia l formu l a for calculating 



airmass (X) from zenith distance (z). In his review, lYoimd (|1974I ^ discussed 



airmass calculations and advocated the use of a simpler formula, 

X = sec z[l - 0.0012(sec2 z - 1)], 
which is valid up to secz = 4 and good to better than 1%. lYound (^O^) also 



warns that airmasses greater than 4 should be avoided because of large random 
and systematic errors. 

Airmass can change significantly during a long exposure. In such cases we 
can use effective airmass (Stetson 1989), 

-^eff = 7 [AsTART + 4AmidDLE + Aend] • 
D 

This is valid provided that extinction is not varying significantly as during the 
exposure. An alternative is to take multiple shorter exposures. 

3.3. Seeing 

In slit-less spectroscopy resolution depends on seeing. We can characterise the 
spectral resolution as follows, 

AA 



^ V ^seeing + ^hist' 

where ^seeing is the size of the stellar seeing disk and 6inst is the fixed instrumental 
resolution for an idealised point-size stellar image. 

Seeing is naturally highly variable, with not only temporal variations (0(t)\ 
but a lso variations with airmass {9 ~ AO-6) and wavelength (6* ~ X'^-^) iMom 
The variation with wavelength can be regarded as fixed, but this resolu- 
tion is different from the instrumental value as determined by arc spectra. 

3.4. Spectral Resolution Issues 

Since the resolution of the spectra is variable, we must allow for this during the 
extinction determination. This is most important where spectrum is varying 
most rapidly, e.g. Balmer lines. Table ^ gives the results that are obtained 
when simple Bouguer law is used to fit to various points within the Balmer 
line during a simulated night of relatively poor seeing. The fit for the Balmer 
core is significantly different from the true value, as can be seen from Fig. |31 
which shows that the depth of the core varies noticeably with airmass. 

The effects of varying resolution become important on nights of poor seeing 
and when seeing is variable. 

3.5. Scintillation 

Scintillation noise is due to motion in the Earth's atmosphere during an exposure 
and sets a lim it on photometric accuracy. This can be estimated using the 
iDravineJ ((199^ formula. 



a = OmD-^/^X'-^'exp ( A ) (2T)-^/^, 
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Table 1. Extinction determined from Bouguer Law fits at various wave- 
lengths within the H/3 Balmer lines compared to the actual value. 
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Figure 3. The effect of resolution varying with airmass on the shape of H/3. 
The core is noticeably affected, and note the more subtle effect in the wings. 

where a is the scintillation noise in magnitudes, D is the diameter of the telescope 
(cm), h altitude of the telescope (m), H scale height of the atmosphere (~8000m) 
and T is the exposure time in seconds. 

The normal minimum exposure time for ASTRA is 10 seconds in order to 
minimise scintillation noise. A 10-second exposure gives a = 0.001 mag. at 
the zenith and 0.007 mag. for X = 3. Longer exposures give lower levels of 
scintillation noise. 

4. Extinction Determination 

Extinction determination will be performed without the use of any assumed 

stellar true fluxes. The extinction stars are tested to be constant, with any 
variability <C 0.01 mag. Any extinction stars found to be significantly variable 
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will be removed from the list of standard stars and extinction re-determined. 
Extinction determination is performed for a single night of observations and its 
reliability assessed. We will investigate the use of multi-night reductions. 

The brightest stars have to be observed through a neutral density filter. We 
will determine and monitor the throughput. The filter transmission function 
enters the extinction deter mination as an extra term in the fitting procedure 
ISterken fc ManfroidI (119921 ). 



Atmospheric extinction is variable. This usually a ppears as a gradual de- 
crease during the night, due to a slow fallout of aerosols ((Yound 119741 ). For a few 
hours, the extinction change can be represented by a linear function with time. 
Observationally, we must observe enough extinction stars to maintain a nearly 
continuous check on the extinction coefficients. Short exposures are required to 
avoi d significant changes in airmass and extinction. 

lYoung Sz Irvind ()l967j ) weighted each observation by 1/ sec z so that resid- 
uals have units of mag. /airmass. Thus residuals against time allows deviations 
from mean extinction to be spotted. 

Initially we assume a perfect night and use a time-independent mul ti-star 
Bouguer Law fitting procedure as described in lSterken Sz ManfroidI ( 1992t ) using 



least-squares techniques ( Pre ss et a l. 1992, Chap. 15). The residuals from the fits 
are assessed using a variety of techniques to look for variations and correlations 
with time and/or airmass. Checks are also made to ensure that none of the 
extinction st ars is variable. As well as e xamining the x^j we will use R-statistics 
as defined bv lBautista &: Steineil ((19931 ). 



^ ^ 1 'EyiVi+iWi 

X;o-2(n-l) Y.Wi/n 

where Wj = — ^ — . These measure the degree of correlation of the residuals for 
adjacent values, which must be strictly monotonic. 

In addition, a 2-d grey-scale image plot of residuals against wavelength and 
time is a very useful diagnostic. For example. Fig. 0] shows a simulation of a 
night with highly variable extinction, including changes in ozone and water, and 
a spell of grey absorption due to a passing cloud band. 

If variable, we can use k{t) as a low-order polynomial in time, adding extra 
terms and using an F-test for the significance of improvement (Bevington 196§, 
Chap. 10). If necessary, we can split the night into smaller blocks, in order to 
avoid parts with poor conditions. 



5. Fitting Telluric Regions 

Non-linear extinction violates Bouguer's Law. Several methods have been used 
in the literature, but we will use synthetic telluric spectra. Non-linear least- 
squares fitting is used to determine oxygen and water columns and their temporal 
variability. Theoretical tra nsmission functions c an be calculated using the HI- 
TRAN molecular line lists if Rothman et al.ll2005^ and a simplified 6-layer model 
of the Earth's atmosphere (Mcholls 198^ . These are calculated with varying 
water and oxygen columns and tabulated at high resolution (0.5A steps) for 
rapid interpolation. 
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Figure 4. Simulation of a night with variable extinction. The residuals from 
the extinction determination reveal obvious variations in ozone and water 
vapour, as well as a spell of grey extinction due to a simulated cloud band 
between frames 140 - 210. The image is ordered by frame number as a proxy 
for time. 



At the resolution of the ASTRA instrument the individual telluric lines are 
not resolved. Thus, flux through each resolution element will depend on the 
nature of the intrinsic telluric lines and any stellar features within that region. 
Strong absorption (or emission) features may be hidden by the lines or 'visible' 
bet ween them. We c annot a priori know which case is occurring in each band 
fe.g. IStevensT^ll 9941 V We also have to beware of stellar lines (dis)appearing in 
and out of telluric lines due to radial velocity variations. The Earth's orbital 
velocity (30 km s~ ) gives a 0.5A shift at 5000A. We wiU assess the level of this 
uncertainty using synthetic stellar spectra. 



6. Flux Calibration 

Spectrophotometric observations must be placed on an absolute flux scale by 
reference to stars with known values of true flux at top of Earth's atmosphere. 
Standard stars have been calibrated against terrestrial sources of known prop- 
erties. Unfortunately, very few stars have been calibrated directly. The pri- 
mary standard is a Lyr ( Vega). The standard referenc es for absolute fluxes are 
iHaves fc LathamI (|l975l V iTiig et al.l (|l977l 'l and Haves (1985). The accuracy of 
these and other calibrations was assessed hv iMeeessieri ((1995|). 
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Most calibrations are at relatively low resolution (typically 50-lOOA). Cor- 
rections for bandpass effec ts need to be considered, which may require the use 
of high-resolution spectra. IColina eiTal . (1996J presented a calibration for Vega, 
for use with HST data, which used Kurucz ATLAS9 fluxes to generate a higher- 
resolution. A m ore-recent absolute calibration based STIS observations was de- 
termined b yiBohlin &: Gilliland (,2004 1 . A discussion of HST standards is given 
bv lBohlinI (^2007^ . 

The ASTRA spectrophotometry will be placed on the Vega absolute flux 
scale. This will require the use of a neutral density filter, whose throughput will 
be measured and periodically checked. Whenever, we use the neutral density 
filter we will observe a suitable star with and without the filter to determine the 
throughput. 

Our quality control procedures will assess stellar (micro) variability and the 
long term variability of extinction stars to ensure they are constant. In addition, 
we will monitor the instrumental and telescope throughput for variations. These 
tests will ensure flux consistency within secondary standards and of the fluxes of 
all target stars. Re-calibration of whole data set be will performed as necessary. 
For standards we will use only the best observations at the highest resolution. 
For other observations, we will ensure spectral resolution effects are allowed for. 

Flux calibrated spectra for all observations will be produced. The individual 
observations of non-standard stars will be co-added, if non-varying, in order to 
increase the signal to noise. Variable stars will not be co-added, except if the 
variations are insignificant between observations. 

The reduction procedures will store intermediate results, such as extracted 
spectra in counts, instrumental counts/second above Earth's Atmosphere, de- 
tails of the nightly extinction models, and individual flux calibrated spectra. 

The final fluxes will be made available as FITS files including headers and 
processing history. For each spectral point, we will give wavelength, bin size, 
instrumental (counts/s), flux relative to Vega, absolute fluxes (photons/s/nm), 
internal errors, external errors, and quality flags. 



7. Conclusions 

The reduction of ASTRA spectrophotometry requires careful extraction and cal- 
ibration in order to achieve stellar flux measurements with internal (star-to-star) 
precision better than 1% for stars brighter than ~10th magnitude. In the regions 
heavily affected by telluric lines we will not necessarily reach the 1% precision 
level. For all observations, including faint spectrophotometric standards, we will 
use full error propagation to include both internal and external errors. 

The available absolute calibrations are accurate to typically 1-2%. Ulti- 
mately this uncertainty will limit the accuracy of the final fluxes of other stars. 
However, the internal precision will be significantly higher, and should more- 
accurate absolute calibrations become available the fluxes can be re-calibrated 
to higher accuracy. 
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